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•  Examples  of  micro  propulsion  ablative  devices 

•  Fundamentals  of  ablation 
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MichiganEngineering 


Istanbul,  Turkey,  June  2004 


2 


Micro-Pulsed  Plasma  Thruster 

(AFRL) 
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Teflon  Propellant 


Plasma 


1 

^  Current 

Electrodes 

PULSER 

1 

Triggerless  ignition 


MicroPPTs 


100  grams 

Energy  1-10  J 

Thrust- to-Power  1-10  juN/W 

Simple  Engineering 

Ablation  rate: 

10  cm  in  100  hour,  1  Hz 


FalconSat 

US  Air  Force  Academy,  2005 
<100-kg  class 
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Micro-PPT  Technology 
Development  (JHU-APL) 
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•APL  exploring  fundamental  effects 
associated  with  device  scaling 

•Examining  influence  of  energy 
deposition  and  electrode  geometry  upon 
electro-mechanical  response 

•Developing  novel  techniques  for  micro- 
PPT  device  characterization 
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Micro-Vacuum  Arc  Thruster 

(AASC) 
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Anode 


Magnetic 

coil 


Inductive  energy  storage  PPU 
(efficiency  of  the  PPU  >90%) 
Low  mass  (<300g) 

Typical  current  -100  A 
Voltages  of -25-30  V. 

Total  efficiency  -10% 

1  -1000-3000  s 
T/P-10  pN/W 
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Micro-Laser  Plasma  Thruster 
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•micro  Laser-ablation  Plasma  Thruster,  p  -LPT  (Photonics  Associates), 
-micro  chip  Laser-ablation  Plasma  Thruster,  (Lincoln  Lab). 


"Reflection  Mode” 

Pro:  2X  higher  Isp,  Cm 
Con:  Optics  Require  Shield 
Awkward  Geometry 


1W  Diode 


Ablatant 


\ 


Transparent  Substrate  /  \  e.g,,  acetate  film 


“Transmission  Mode” 

(Analogous  to  LIFT) 
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Ablation  Fundamentals 
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Layer  structure  near 
ablated  surface 


Istanbul,  Turkey,  June  2004 


7 


Kinetic  model  of 
the  Knudsen  Layer  ( 1 ) 
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Analytical  and  particle  (dsmc) 
approaches: 


JVxf(V)  dV  =  const 


f(x,V)  =  ^(x)f1(V)  +  (1  -W2(V) 
where  ^(x=0)=1  and  ^(L)=0  with  x=0 
[Mott-Smith,  1951] 


JVx2f(V)dV  =  const 
jVxV2  f(V)dV  =  const 


(mass) 

(momentum) 

(energy) 


f,(V)  =  n0P3/2exp(-V2)  Vx>0 

f-i  (V)  =  Sf2(V)  Vx<0 

f2(V)  =  n1p3/2exp(  -(v-U)2) 
[Anisimov,  1968] 


— ■ —  analytic 
— DSMC(U) 
—a—  DSMC  (10A.) 
DSMC(IOOA) 


0.5 

U1/(2kT1/m)° 
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Kinetic  model  of  the 
Knudsen  Layer  (2) 


Particle  distribution  function.  DSMC  Parameters  at  the  Knudsen  layer  edge 


U1/(2kT1/m)05 
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Hydrodynamic  Layer 
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d(nV)  =  0 

.  .  dV  d(nkT) 

M(nV) - = - - - -  +  jxB 

dx  dx 

nt  T2n2  |  1  ju(jd)2 
MV?  _  2  27;  4  kTj 

2kT.  ~  3  n.2 

- -  - 

2  «2 


Vj  depends  on  the  specifics  of 
acceleration  (n2,  j) 


Example:  electrothermal  thruster 
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Velocity  at  the  Knudsen  layer  edge 


Velocity  at  the  Knudsen 
Layer  Edge 
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Discharge  energy  (J) 

(Current  density) 

Electromagnetic  acceleration  leads  to 
transition  to  vacuum  evaporation  regime 


Mikellides  et  al 


Electrothermal  PPT 


Keidar  et  al  2000 

In  addition  to  acceleration  mechanism 
PPTs  can  be  classified  by  ablation  mode 
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End-to-end  simulation 
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1/4”  DIA  (J-PPT,  6  J 


Comparison  with 
experiment  ( 1 ) 
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1 3.2  -  10 
j  1.0  -  3.2 
0.32  -  1.0 
0.10  -  0.32 
0.03  -  0.10 
0.010  -  0.03 
0.003  -  0.010 
I  IE-3  -  0.003 


a  z,  mm 

6  mm  beam  width 

Herriot  Cell  interferometer 
Antonsen,  Burton  &  Spanjers,  IEPC-01-157 
3  mm  Antonsen,  MS  Thesis,  UIUC  2001 
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Electron  density  (xIO' 


Comparison  with 
experiment  (2) 


MichiganEngineering 


1/4”  DIA  p-PPT,  6  J 
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Cathode 


Propellant  recession  (1) 
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■ 


Anode 


2D  PIC-DSMC  model  & 
magnetic  transport 

-  time-dependent 
boundary  conditions 
plasma  layer  model 

-  magnetic  field  & 
current  distribution 
(energy  balance  &  ion 
dynamics) 

-  collisions  (elastic  & 
non-elastic) 
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Propellant  recession  (2) 
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1/4”  DIA,  6  J 
6  hours,  1Hz 
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Ablation  pattern: 
Charring 


MichiganEngineering 


Ablation  rate  calculations  are 
based  on  plasma  layer 
model  and  ablation  theory 


AFRL  Experiment: 

Tests  agree  with  model  - 
Cold  thrusters  char  easier 
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Discharge  non-uniformity 
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High  speed  camera 
visible  emission 


Arc  spoking  increases  with  energy 


2.30  J  5.55  J  6.73  J 
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Radial  distance  (m) 


Current  constriction  modeling 
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min  ii 


Surface  temperature  (K)  Surface  temperature  (K) 


1200-1 


Surface  temperature  ( 1 ) 
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Experiment  ill  the  Spot 


13AUG  Temperature  at  30  us 
1 0  J  test 


Shot  number 
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Micro-Laser  Plasma  Thruster 
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Plasma  Jet 
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Energy  equations: 

VdTp/ck  =  Q1B  -  Qei  -  Qx 

Dominated  by  inverse-Bremsstrahlung 

Q w  —  ^ w  Li  cxp(— CCjgX) 

J 

ccIB  =  1. 37  xlO~35  A3  Ne2Te  2 

Plasma  properties: 

-  composition  (C,H,C+,H+,  qj 

-  assuming  LTE 
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Radial  Distance  (m) 


Plume  Simulation  (PIC-DSMC) 
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Comparison  with  experiment  ( 1 ) 
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Experimental  set  up 
Witness  plate  deposition 
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Comparison  with  experiment  (2) 


Ql-lQ-SOsi 


SEM  analysis  shows  that  deposition  material  is  carbon 
Thus  carbon  flux  is  compared  with  deposition  profiles 
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Comparison  with  experiment  (3) 

2.5  W;  P=10-2  torr  8  W;  P=6.5xlQ-5  torr 
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Background  pressure  and  operational  conditions 
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Micro- Vacuum  Arc  Thruster 


MHD  free  boundary  model 


mi  (V.-VJVj  =  -k(Z.  T+T^VMn)  +  j  x  B/;; 
j  =  a{E+(/tre  !e)'V\n{n)  -  jxB /(en)  +  (V,xB)| 
V-(V,  n)  =  0 
V-j  =  0 
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Density  Density  Density 


Micro- Vacuum  Arc  Thruster 

Plume 
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Summary  ( 1 ) 
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-  Various  microthruster 
technologies  based  on 
ablative  mechanism  were 
developed 
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Summary  (2) 


-  Self-consistent  modeling  approach  for  ablative 
micro-thrusters  was  formulated  based  on  a 
kinetic  ablation  model  and  particle  plume 
simulation. 

-  Most  extensive  validation  of  the  modeling 
approach  was  performed  for  micro-pulsed 
plasma  thruster.  Plasma  density,  surface 
temperature,  ablation  rate,  ablation  profde  were 
compared  with  experiment.  Optimization  criteria 
were  formulated  for  some  devices,  such  as 
microPPT. 
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Future  needs 
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•  Development  of  more  flexible  technology  (variable  Isp,  variable  thrust) 

-  pLPT  pulse  duration 

-  pVAT  pulse  duration,  material 

•  Contamination  issues 

-  Study  is  needed  (pVAT) 

•  Lifetime  issues  (propellant  recession) 

-  pVAT,  pPPT 

•  Hybrid  thrusters 

-  pVAT/pPPT;  pLPT/pPPT 

•  Modeling:  further  characterization  of  thrusters  (pVAT),  plumes.  Effect 
of  the  magnetic  field 
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